Partly Occupied Wannier Functions 
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We introduce a scheme for constructing partly occupied, maximally localized Wannier func- 
tions (WFs) for both molecular and periodic systems. Compared to the traditional occupied WFs 
the partly occupied WFs posses improved symmetry and localization properties achieved through 
a bonding-antibonding closing procedure. We demonstrate the equivalence between bonding- 
antibonding closure and the minimization of the average spread of the WFs in the case of a benzene 
molecule and a linear chain of Pt atoms. The general applicability of the method is demonstrated 
through the calculation of WFs for a metallic system with an impurity: a Pt wire with a hydrogen 
molecular bridge. 

PACS numbers: 71.15.Ap, 31.15.Ew, 31.15.Rh 



The idea of representing the electronic ground state of 
a solid-state system in terms of localized orbitals dates 
back to Wannier [l| who introduced a canonical transfor- 
mation between the Bloch states of an insulating crystal 
and a set of localized functions. These localized Wan- 
nier functions (WFs) provide a formal justification of the 
widely used tight-binding 0] and Hubbard models Q • A 
scheme for evaluation of maximally localized WFs for 
composite bands has been developed by Marzari and 
Vanderbilt Q, and the WFs have turned out to be use- 
ful in a number of contexts including electron transport 
calculations ||, linear-scaling electronic structure meth- 
ods and in the theory of electric polarization |3, la] ■ 

In the context of molecular systems the analogous of 
Wannier functions for finite systems has been studied un- 
der the name "localized molecular orbitals" 0, flfl ITU 
E 13 • These are traditionally defined by an ap- 
propriate unitary transformation of the occupied single- 
particle eigenstates and represent a strong analysis tool 
as they provide insight into chemical properties such as 
bond types, coordination, and electron lone pairs. In the 
following we shall for simplicity use the term WF to cover 
also localized molecular orbitals. 

In this letter we demonstrate that in constructing lo- 
calized WFs it is in some cases advantageous to con- 
sider not only the space of occupied states but to in- 
clude also selected unoccupied states through a bonding- 
antibonding closing procedure. The changes achieved by 
the bonding-antibonding closing can be viewed directly 
in Fig. ^ where calculated WFs with and without the in- 
clusion of additional unoccupied states are shown in the 
case of a benzene molecule. In the upper panel only the 
occupied states are used and some of the resulting WFs 
clearly mix a- and 7r-character. By including three se- 
lected unoccupied states the resulting partially occupied 
WFs become more localized and the a- and 7r-systems 
are now separated. In addition the rotational symmetry 
of the molecule is more clearly reflected by the partly 



occupied WFs. 

The problem solved by the bonding-antibonding clos- 
ing is the following: if we consider two well-localized 
atomic orbitals on neighboring atoms in a molecule and 
allow the two states to hybridize, a bonding and an an- 
tibonding combination will result - combinations which 
may be less localized than the individual atomic orbitals. 
To regain the localized atomic orbitals from the molec- 
ular orbitals we need both the bonding and the anti- 
bonding combination independent of their occupation. 
In some cases the antibonding state may have hybridized 
further with other states and the state which "matches" 
the bonding state will be a linear combination of eigen- 
states. To apply the bonding-antibonding closing in prac- 
tice a systematic scheme to identify the relevant unoccu- 
pied states is clearly required. We show that this can be 
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FIG. 1: Occupied and partly occupied Wannier functions 
for benzene CeH6. The traditional (occupied) WFs are seen 
to mix the a and % orbitals. Inclusion of the three anti- 
bonding 7r orbitals leads to the partly occupied WFs which 
are more localized and separate the orbitals with different 
mirror symmetry. 
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achieved by optimizing the localization of the resulting 
WFs. 

The lack of localization of WFs obtained if only oc- 
cupied states are considered is well-known in the case of 
extended metallic systems 0, 0. A proper localiza- 
tion of the WFs requires in this case completion of the 
partly occupied valence bands (the analogue of bonding- 
antibonding closing for a periodic system) which may be 
difficult due to crossings and mixing of the valence bands 
with unwanted higher-lying bands. This problem was ad- 
dressed by Souza, Marzari, and Vanderbilt who con- 
structed a disentangling procedure which makes it possi- 
ble to follow a given band by minimizing the change in 
character of the Bloch states across the Brillouin zone. 
In this way the relevant unoccupied states could be iden- 
tified. The approach we suggest here makes no reference 
to the wave vector, it is simple to apply and valid for 
both isolated and periodic systems. 

We note that different techniques for constructing lo- 
calized orbitals have been developed in connection with 
linear scaling electronic structure methods [||. In that 
context the orbitals are typically constructed during an 
electronic structure calculation and are required to vanish 
outside predefined localization regions. In contrast, the 
scheme we propose is an exact transformation of already 
determined eigenstates to obtain maximal localization. 

Our starting point is a set of single-particle eigenstates 
{ipn} obtained for example from a Density Functional 
Theory (DFT) 17] calculation. We denote the total num- 
ber of eigenstates included in the DFT calculation by N 
and the number of occupied eigenstates by M . Following 
Marzari and Vanderbilt [4j we measure the spread of the 
WFs {w n } by the sum of second moments 

S = ^2((w n \r 2 \w n ) - (w n \r\w n ) 2 ). (1) 
n 

The WFs are expanded in terms of the M occupied 
eigenstates and L extra degrees of freedom (EDF) 
from the (N — M)-dimensional space of unoccupied eigen- 
states. This will lead to a total of M + L WFs and we 
thus require N > M + L. The expansion takes the form 

M L 

W n = 2J U mn 1p m + 2J UM+l,n<Pu ( 2 ) 
rn=l 1=1 

where the unoccupied orbitals are written as 

JV-M 

4>l = X! c ml'>pM+m- (3) 
m— 1 

The columns of the matrix c are taken to be orthonormal 
and represent the coordinates of the EDF with respect 
to the unoccupied eigenstates. The matrix U is unitary 
and represents a rotation in the space of the functions 
{ipi, ■ ■ ■ ,ipM, 01 1 • • • j <t>t}- Through the expansions 



and , S becomes a function of E/y and Cy , which must 
be minimized under the constraint that U is unitary and 
the columns of c are orthonormal. 

To keep technicalities at a minimum the formulation 
given above has been specialized to systems which are 
either isolated or described in a periodically repeated su- 
percell in the T-point approximation. A detailed account 
of the method including the k-point formulation will be 
given elsewhere [l8| . Here we focus on the case of a cubic 
supercell of length a in the T-point approximation. Ex- 
tension of the theory to supercells of general symmetry 
is straightforward along the lines of Refs. andfilil 

In the limit of a large supercell, minimizing S is equiv- 
alent to maximizing the functional |13l |20j 

fi = J2(\ X ™\ 2 + \Ynn\ 2 + \Z nn \ 2 ), (4) 
n 

where the matrix X is defined as X nm — 
(w n \e~ l ^/ a ^ x \w m ), with similar definitions for Y nm and 
Z„ m , The unitary matrix at iteration n is written as 
{/(«) = [/( n_1 )exp(— A), where A is an anti-hermitian 
matrix. By straightforward differentiation we obtain 
expressions for the gradients dCl/dAij and dCl/dc*^ 
giving the steepest uphill direction of Cl upon variations 
in A and c. The orthonormality constraint on c is 
invoked through a set of Lagrange multipliers. Any 
gradient-based optimization scheme can now be applied 
to maximize Cl. 

Given N and L, the algorithm introduced above pro- 
duces the M + L most localized WFs. We shall assume 
that N has been fixed at a value large enough to include 
all anti-bonding states relevant for the localization. The 
problem is then to determine a good value for L. It 
seems as a natural strategy to choose an L that gives a 
high localization per orbital. To quantify this statement 
we define the average localization per orbital as 

(fi) = (5) 
M + L K ' 

where M + L is the total number of WFs. From the 
definition of £1 it is clear that < (O) < 3. Fixing L 
on the basis of (Q) represents a quite general criterion 
which can be applied to any system. However, the local- 
ization procedure must be carried out for several L which 
might be a tedious task depending on the size of the sys- 
tem. Formally the global maximum of (CI) is attained in 
the limit where both ./V and L tend to infinity in which 
case completely localized delta functions can be realized. 
However, we have found that for practical values of N 
where very high energy states are not considered, (CI) 
goes through a local maximum as L is increased beyond 
0. For many systems it is possible to determine a value 
for L based on symmetry arguments, chemical intuition, 
or a closed band condition. In particular the latter strat- 
egy is natural for crystalline systems where the closed 
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FIG. 2: Average localization per Wannier function as a func- 
tion of the number of included unoccupied orbitals for the 
three systems studied. For the benzene molecule the maxi- 
mum appears at L — 3 corresponding to the inclusion of the 
three antibonding n orbitals. For the Pt chain the maximum 
at L = 6 corresponds to completion of the valence bands. 
(The values of (fl) do not compare from one system to an- 
other due to dependence on the size of the supercell.) 



band condition is satisfied by choosing L equal to the 
number of unoccupied states belonging to all partly oc- 
cupied bands. As we shall see in the following examples 
the two criteria for determining L lead to similar results. 

All the following calculations are based on the Kohn- 
Sham eigenstates as defined within DFT We use a 
plane wave based pseudopotential code [2l| with an en- 
ergy cutoff of 25 Ry for the plane wave expansion and de- 
scribe the ion cores with ultrasoft pseudopotentials ■ 

We have constructed the maximally localized WFs of 
an isolated benzene molecule in a cubic supercell of length 
a = 16 A. If we consider only the occupied eigenstates of 
benzene, corresponding to L = in the method, we find 
the 15 occupied WFs shown in the upper panel of Fig.Q 
Nine of these are cr-bonds centered at each C-H and ev- 
ery second C-C bond, and the last six orbitals are double 
bonds with a mixed a- and 7r-character centered at the 
three remaining C-C bonds. This is the classical picture 
of the chemical bonds in benzene. In this picture the re- 
sulting orbitals apparently break the six-fold rotational 
symmetry of the molecule, i.e. the bonding orbitals be- 
tween different C-C pairs differ. Moreover, the planar 
structure of the molecule suggests a separation of a- and 
7r-orbitals contrary to C-C double bonds. The occupied 
subspace spanned by all the WFs of course posseses the 
full symmetry. 

The lack of symmetry in the chemical picture of ben- 
zene is a result of ignoring the three anti-bonding ir- 
orbitals of the C ring in the construction of the WFs. 
In Fig. |21 (triangles) we plot the average localization of 
the WFs for various L values. The maximum occurs 
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FIG. 3: The norm of the projection of eigenstates onto the 
subspace spanned by the partly occupied WFs for the three 
systems studied. The eigenstates are listed in increasing order 
according to their energy. States below the dashed line are 
occupied. It is seen that in the case of the benzene molecule 
and the Pt chain the algorithm selects particular unoccupied 
states necessary for the bonding-antibonding closing. In the 
case of the hydrogen molecule in the Pt chain one of the WFs 
has weight on an unoccupied resonance distributed over sev- 
eral eigenstates. 



for L = 3. To analyze the composition of the EDF we 
have projected the molecular eigenstates onto the WFs 
generated with L = 3, see upper panel of Fig- EI By m_ 
spection we have verified that the EDF selected by the 
algorithm are exactly the three anti-bonding 7r-orbitals 
of the C ring. Maximizing (f2) is therefore equivalent to 
performing a bonding-antibonding closure. The 18 partly 
occupied WFs obtained with L = 3 consist of cr-bonds 
at all C-H and C-C pairs and a p-orbital perpendicular 
to the molecular plane centered at each C atom (lower 
panel of Fig. Q}. Consequently, the inclusion of 3 EDF 
restores the rotational symmetry and separates a- and 
7r-orbitals. 

We have performed a similar analysis for a linear chain 
of six Pt atoms with periodic boundary conditions (Ptg), 
see Fig. Hfla). Such monatomic chains have been pro- 
duced and studied in STM experiments |2^| and by the 
mechanically controlled break junction technique [24| . 

The valence bands of an infinite Pt wire consist of five 
d- and one s-band. For each band only six Bloch states 
comply with the boundary conditions in Ptg and conse- 
quently it takes a total of 36 states to close the bands. 
Since there are 30 occupied states, the closed band con- 
dition can thus be fulfilled with L = 6. From Fig. |3 
(circles) we see that the closed band condition coincides 
with the maximum of (J7). In the middle panel of Fig.|3 
we show the projection of eigenstates onto the Wannier 
functions corresponding to L = 6. By comparing with 
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FIG. 4: Partly occupied Wannier functions for a periodically 
repeated, linear string of 6 Pt atoms without and with an 
inserted hydrogen bridge. The WFs have been generated with 
L = 6 (Pt 6 ) and L = 7 (Pt 6 -H 2 ). The WFs shown for the 
Pt6-fi2 system correspond to H2 and H-Pt <r-bonds. 



be expanded in terms of the obtained WFs. The inclusion 
of additional unoccupied states in the definition of the 
WFs improves the localization and simplifies the result- 
ing chemical pictures. The method can easily be modified 
so that not the occupied states but for example the states 
in an energy window around the Fermi level are exactly 
represented. Extra high- and low-energy states can then 
be included to construct effective model Hamiltonians for 
the states around the Fermi level. 

We acknowledge discussions with S0ren L. Frederiksen 
and support from the Danish Center for Scientific Com- 
puting through Grant No. HDW-1101-05. 



the band diagram of the Pt wire, we find a perfect match 
between the EDF and the unoccupied Bloch states of 
the valence bands. Again this shows that closing the 
bands is equivalent to minimizing the average spread. For 
L = 6 we find five atom-centered d-orbitals, of which one 
is shown, and one cr-orbital centered between each Pt- 
Pt pair. The transformation to WFs allows for explicit 
extraction of tight-binding parameters. For example, the 
doubly degenerate (5-band of the Pt chain is characterized 
by the hopping parameters (in eV) t\ = 1.08, i2 = 0.13, 
t 3 = 0.03. 

When N is increased from 50 to 90 the localization gen- 
erally improves due to the larger space available for the 
EDF, see Fig. [21 (open circles). If we increase N further 
this trend would continue and eventually the maximum 
would start shifting towards larger L values. In the limit 
N — > 00 (f2) would increase monotonically as a function 
of L. This shift of the maximum, however, occurs only 
when N becomes very large, i.e. when high-energetic and 
chemically irrelevant states enter the localization space. 

As a final example we analyze the effect on the WFs 
when a molecular hydrogen bridge is inserted in the Ptg 
system, see Fig.|3f2a). A molecular hydrogen bridge be- 
tween Pt contacts has recently been realized and studied 
experimentally j25T ] . The average localization has a max- 
imum for L — 6. The maximum is, however, not very 
pronounced and its position depends on the Pt-H and H- 
H bond lengths. Consequently we have studied the WFs 
for L = 5, 6, 7. In all three cases we recover the five d- 
orbitals per Pt found for Pt6- The er-orbitals are slightly 
distorted and the one in the bridge region is replaced by: 
an H2 bonding orbital (L — 5), two s-orbitals located 
symmetrically on each H atom (L = 6), an H2 bonding 
orbital and two Pt-H bonds with a symmetry (L = 7). 
The latter situation is shown in Fig. Ef2b). As can be 
seen from the lower panel of Fig. [3] this system provides 
an example where one of the EDF is a superposition of 
eigenstates. 

To conclude we have proposed a method for finding 
maximally localized, partly occupied WFs. The method 
ensures by construction that any occupied eigenstate can 
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